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SUMMARY 


This paper describes the development and usage of the 
Thermal Energy Flow Simulation (TEFS) computer program. The 
program's function is to simulate transient heat conduction in 
composite solids and accurately predict interfacial temperatures. 

Development of this program was undertaken in order to 
analytically test the thermal characteristics of a proposed 
thermal protection tile system for the Space Shuttle Orbiter. 
This tile configuration (referred to as a Shell Tile*) consists 
of an outer metallic shell which is filled with multi-layered 
flexible insulation. The objective was to determine individual 
transient temperatures between layers of dissimilar insulation 
materials and to accurately predict the peak temperature of the 
underlying thermally protected aluminum structure. A major 
assumption was that (along with each material layer's thermal 
properties) the only known variable is the temperature of the 
outer surface. 

The thermal protection tile system was heat tested and a 
time history of the interfacial temperatures was recorded. Using 
the TEFS program, a computer simulation was conducted to 
analytically determine the temperature distribution. The 
computer simulation provided a peak aluminum temperature that was 
within 4% of the actual. 


♦United States patent number 4,456,208 dated June 26, 1984 

entitled "Shell Tile Thermal Protection System". 
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INTRODUCTION 


A new thermal protection system, referred to as a Shell 
Tile, has been proposed for possible application to future entry 
vehicles. In order to optimize the tile thickness via an 
analytical model, the Thermal Energy Flow Simulation (TEFS) 
computer program was developed. The unique feature of this 
program is its ability to accurately simulate the thermal 
characteristics of a composite solid when the only data available 
are the top surface time dependent temperature and the materials' 
thermal properties. 

The basis for the TEFS program is the THERMAL RESPONSE 
FACTOR (RESFAC) program used in NASA'S. ENERGY COST ANALYSIS 
PROGRAM (NECAP) . In both TEFS and RESFAC, the equation of heat 
conduction for unsteady state/linear flow is solved using the 
Laplace transformation method. This method was used due to its: 

1) Direct application to composite solids. 

2) Easy input requirements. 

3) Fast computer execution time. 

A detailed analysis of the method can be found in references 1 
and 2 . 

Note that the following assumptions govern the use of the 
Laplace transformation method as used in TEFS: 

1) Heat transfer is by conduction in one dimension. 

2) No contact resistance exists between f mater i al layers. 

3) The composite solid is initially at a uniform 
temperature. 

4) Values of temperature dependent thermal properties are 
determined by the temperature at the top of the 
respective material layer. 

5) The last layer of the composite solid is either an air 
or vacuum gap which is held constant at the initial 
temperature . 

In order to obtain an accurate simulation, the following 
criteria should be observed: 

1) Thick materials should be subdivided. 

2) The last material of the model being analyzed should be 
an air/vacuum gap whose only thermal property is a 
constant resistance. 

3) Any material which has negligible thermal mass should 
have resistance as its only thermal property. 

4) The time step selected for the simulation should 
accurately reflect the supply temperature's duration. 

5) Temperature dependent thermal properties should be input 
such that the temperature range being simulated is 
accurately represented. 



PROGRAM INPUT/OUTPUT 


Two input files are required for the TEFS program. The 
first file, accessed under file name MODEL, contains the tempera- 
ture dependent thermal properties of each material layer compris- 
ing the composite solid. File two, accessed under file name 
SUPPLY, contains values defining the thermal simulation. Note 
that data in both files are free formatted. 

The contents of file MODEL are: 


Item 


Description 


N, IU 


N: integer number of material 

layers; maximum of 10 
IU: integer units flag 

0 for English, 

1 for Metric 




,1 


N 


Integer number of temperature 
variations in layers 
1 ( top) , 2 ,... ,N ( bottom) ; maximum 
of 20 for each layer; note that 
a layer having a constant 
resistance must be entered as a 
-1 


NOTE: Repeat the following for all material layers, proceeding 

from top layer to bottom 


LABEL. 


Type of material used in the 
ath layer; alphanumeric data 
with a maximum of 80 characters 


NOTE: The following input is required for each temperature 

variation in the ath layer 


W,K,D,S,T 


If layer has thermal mass, 
enter real values for: 
thickness W, density D, 
conductivity K, specific 
heat S, temperature T 


or 


R 


Layer has no thermal mass, 
enter real. value for: 
constant resistance R 


The contents of file SUPPLY are: 


Item 


Description 


DT Real value of time increment 

between supply temperature 
pulses 



Item 


Description 


NP Integer value of number of 

supply temperature pulses; 
maximum of 100 


TS ^ , TS2 / • • • » TSjgp 

The output file from TEFS 
following : 


NP real values of supply 
temperature pulses from TS 
initial through TS end 

is named SIMUL and contains the 


1) Echo of input file MODEL 

2) Echo of input file SUPPLY 

3) Tabular presentation of simulated temperatures at the 
top of each material layer for the number (NP) of time 
steps (DT) . 
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PROGRAM UNITS 


Depending on the value of the units flag IU, the appropriate 
units for the variables in the input files MODEL and SUPPLY and 
the output file SIMUL are: 


Item 

W 

D 

K 

S 

T 

R 


English 

ft 

lb/ft 3 
Btu/h-f t-F 
Btu/lb-F 
F 

h-f t 2 -F/Btu 


Metric 

m 

kg/m 3 

W/m-K 

W-s/kg-K 

K 

m 2 -I</W 


DT 


s 


s 


TS 


F 


K 
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ANALYTICAL VS RECORDED DATA 


A new design for a thermal protection tile has been 
investigated as a possible alternative to the current space 
shuttle orbiter tile system or as a heat shield for future 
aerospace transports. This new design differs significantly from 
« the currently used tiles. Rather then use rigidized ceramic 

tiles, the proposed design consists of individual metallic shells 
which are filled with layers of low density flexible 

* insulation. Figure 1 contains a description of this metallic 

shell concept. 

Heating experiments have been conducted on this tile design 
in order to obtain temperature profiles of the materials' 
interface regions. Using program TEFS, an analytical solution 

was obtained. The TEFS input files MODEL and SUPPLY are listed 
in TABLES I and II respectively; output file SIMUL is presented 
in TABLE III. A graphics comparison between recorded and 
analytical data is shown in Figures 2 through 7. 

The computer simulation, run interactively on a Control Data 
Corporation CYBER 160/170 series computer, required 40 CPU 
seconds and 73000 (octal) 60-bit words of memory to execute. 

Probable causes of the discrepencies between recorded data 
and the analytical solution are: 

1) Computer simulation is one dimensional heat conduction. 

2) Thermal properties are ideal. 

3) There is an initial temperature gradient within the 
model. 

4) Uneven heating rate exists due to staggered pattern of 
heating elements. 

5) Thermocouple inaccuracies, especially at low temper- 
atures. 
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15 x 15 cm 
Upper shell (ceramic coatings) 
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Self-locking 
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fasteners 


Bond/insulator 


Figure 1. Metallic Shell Tile concept. Numbers refer to 
thermcouple locations and also comparisons with 
simulated temperatures. 
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TABLE I 


INPUT FILE MODEL 


4 


7 0 

17 17 17 17 84-1 


TOP 

RENE 


0. 00042 

5. 1993 

512.0 

0.00042 

5.8925 

512.0 

0. 00042 

6.5858 

512.0 

0.00042 

7.279 

512 . 0 

0.00042 

7.9723 

512.0 

0. 00042 

8 . 6655 

512 . 0 

0. 00042 

9. 1709 

512.0 

0.00042 

9.6763 

512 . 0 

0.00042 

10. 182 

512.0 

0.00042 

10.687 

512.0 

0.00042 

11.193 

512.0 

0.00042 

1 1 . 698 

512.0 

0. 00042 

12.203 

512.0 

0.00042 

12.709 

512 . 0 

0. 00042 

13.214 

512 . 0 

0.00042 

13.72 

512.0 

0.00042 

14.225 

512.0 

TOP 

SHELL QUARTZ INSUL 

0.0625 

0.0108 

3.0 

0.0625 

0.0117 

3.0 

0. 0625 

0.0154 

3 . 0 

0 .. 0625 

0.019 

3 . 0 

0.0625 

0 . 0227 

3.0 

0.0625 

0.0263 

3 . 0 

■ 0.0625 

0.03 

3.0 

0.0625 

0. 03371 

3.0 

0.0625 

0. 0383 

3.0 

0 . 0625 

0.0428 

3 . 0 

0.0625 

0.0496 

3 . 0 

0.0625 

0.0564 

3 . 0 , 

0.0625 

0.0632 

3 . 0 

0.0625 

0.07 

3 . 0 

0.0625 

0.0788 

3 . 0 

0 . 0625 

0.0875 

3.0 

0 .. 0625 

0. 0963 

3 . 0 

MIDDLE SHELL C 

2 LJARTZ IN 

0.0625 

0.0108 

3 . 0 

0.0625 

0.01 17 

3.0 

0 .. 0625 

0.0154 

3 . 0 

0.0625 

0.019 

3 „ 8 

0.0625 

0 . 0227 

3 . 0 

0.0625 

0 . 0263 

3.0 


0.0987 

0.0 

100 . 0 

0 . 103 

0.0 

200. 0 

0. 1072 

0.0 

300.0 

0. 1114 

0.0 

400.0 

0. 1157 

0.0 

500 . 0 

0 . 12 

0.0 

600 . 0 

0. 1242 

0.0 

700.0 

0. 1285 

0.0 

800 . 0 

0. 1327 

0.0 

900 . 0 

0 . 137 

0 . 0 

1000.0 

0. 1412 

0.0 

1100. 0 

0. 1455 

0.0 

1200.0 

0. 1492 

0.0 

1300.0 

0. 1521 

0.0 

1400.0 

0. 1549 

0.0 

1500.0 

0. 1577 

0.0 

1600. 0 

0. 1605 

0.0 

1700.0 

0 . 18 

0.0 

100 . 0 

0.2 

0.0 

200.0 

0.215 

0.0 

300 . 0 

0.23 

0 . 0 

400.0 

0.24 

0 . 0 

500.0 

0 . 25 

0.0 

600.0 

0.255 

0.0 

700.0 

0.26 

0.0 

800.0 

0 . 265 

0.0 

900.0 

0 . 27 

0.0 

1000.0 

0.275 

' 0.0 

1100 . 0 

0.28 

0 . 0 

3 200.0 

0.285 

0.0 

1300. 0 

0 . 29 

0.0 

3 400.0 

0.295 

0.0 

1500. 0 

0. 3 

0 . 0 

1600. 0 

0.305 

0 , 0 

1700.0 

0 . 18 

0 . 0 

100 . 0 

0. 2 

0 . 0 

200. 0 

0 „ 2 1 5 

8 . 8 

7.00 . 0 

0 .. 23 

8 „ 8 

4 08 . 8 

0 . 24 

8 . 0 

580 . 0 

0 . 25 

8 . 0 

600 . 0 
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TABLE I. (cont.) 


0.0625 

0.03 

3.0 

0 . 255 

0.0625 

0.03371 

3.0 

0.26 

0.0625 

0.0383 

3. 0 

0.265 

0.0625 

0.0428 

3.0 

0.27 

0.0625 

0.0496 

3.0 

0. 275 

0.0625 

0.0564 

3.0 

0.28 

0. 0625 

0. 0632 

3. 0 

0 . 285 

0.0625 

0.07 

3.0 

0.29 

0.0625 

0. 0788 

3.0 

0.295 

0.0625 

0.0875 

3. 0 

0. 3 

0. 0625 

0.0963 

3.0 

0.305 

BOTTOM SHELL 

QUARTZ IMSIJL 


0.0625 

0. 0108 

3.0 

0. 18 

0.0625 

0.0117 

3.0 

0.2 

0. 0625 

0. 0154 

3.0 

0.215 

0.0625 

0.019 

3.0 

0 . 23 

0. 0625 

0.0227 

3.0 

0.24 

0.0625 

0.0263 

3.0 

0 „ 25 

0.0625 

0. 03 

3.0 

0.255 

0.0625 

0 . 0337 1 

3. 0 

0.26 

0. 0625 

0.0383 

3. 0 

0 . 265 

0.0625 

0. 0428 

3.0 

0. 27 

0 .. 0625 

0 . 04 96 

3 . 0 

0.275 

0.0625 

0. 0564 

3. 0 

0 . 28 

0. 0625 

0.0632 

3. 0 

0 . 285 

0.0625 

0.07 

3.0 

0. 29 

0. 0625 

0.0788 

3.0 

0 . 295 

0.0625 

0. 0875 

3.0 

0.3 

0. 0625 

0. 0963 

3.0 

0.305 

NOMEX 

FELT 

(TIP) 


0.0133 

0.0198 

5. 4 

0.312 

0. 0133 

0.0232 

5.4 

0 . 32 

0.0133 

0 . 0265 

5. 4 

0 . 335 

0. 0133 

0.0304 

5.4 

0. 338 

0.0133 

0.0347 

5. 4 

0. 342 

0.0133 

0.039 

5. 4 

0.345 

0.0133 

0.0434 

5 . 4 

0.35 

0.0133 

0.0478 

5 . 4 

0.36 

ALUMINUM STRUCTURE CONDUCT ION 

0 . 02 

78.0 

175,. 0 

0 . 202 

0. 02 

84.0 

1 75 . 0 

0 . 2 1 

0 . 02 

90 . 0 

175. 0 

0. 223 

0. 02 

97 „ 2 . 

1 75 . 0 

0 232 


AH IM I KILIM RAD I AT I OKI T O VACUUM 


0.0 

700.0 

0.0 

800.0 

0.0 

900. 0 

0. 0 

1000. 0 

0.0 

1100.0 

0. 0 

1200.0 

0.0 

1300.0 

0.0 

1400.0 

0.0 

1500.0 

0.0 

1600.0 

0.0 

1700. 0 

0.0 

100. 0 

0. 0 

200.0 

0.0 

300.0 

0.0 

400.0 

0. 0 

500.0 

0. 0 

600.0, 

0. 0 

700.0 

0.0 

800. 0 

0.0 

900.0 

0.0 

1000.0 

0. 0 

1 1 00 . 0 

0.0 

1 200 . 0 

0.0 

1300. 0 

0. 0 

1 400 , 0 

0. 0 

1500.0 

0.0 

1600. 0 

0.0 

1700.0 

0. 0 

100. 0 

0/ » 0 

200.0 

0. 0 

300.0 

0.0 

403.0 

0. 0 

500.0 

0.0 

600.0 

0.0 

700.. 0 

0.0 

800. 0 

0. 0 

100.0 

0 ■■ 0 

700. M 

0 , 0 

300 . 0 

0 0 

400. (T 


4 6 ,. 0 


8 



TABLE II. INPUT FILE SUPPLY 


60.0 
28 

65.25 186.7 795.1 706.4 1065.8 994.8 1230.0 1478 . 
1637.9 1654.9 1663.6 1647.8 1647.4 1652.6 1658.4 
1477.6 1293.0 1070.0 915.4 809.7 731.9 670.7 620 . 
544.5 514.6 


4 1614 
1655.9 
8 579 . 



TABLE III. OUTPUT FILE SIMUL 



— 

ECKQ OF 

INPUT FILE 

MODEL 




— 

— - ECHO OF 

INPUT FILE 

SUPPLY — 



TIME 


TEMPERATURE A 

T TOP OF LAYER 



SEC 

1 


T 

4 

5 

6 

0 

65. 250 

65.250 

65 . 253 

65.250 

65 . 250 

65. 250 

60 

186.700 

1 86 .68? 

65. 629 

65.250 

65. 250 

65. 250 

120 

795. 100 

795. 037 

72.873 

65.769 

-65 ,. 305 

65. 250 

180 

706., 400 

706. 795 

375. 139 

65.73 0 

65 , 265 

65. 257 

240 

1 065 . 800 

1065. 753 

298 39 1 

68 ., 958 

-65 ,. 3 1 4 

l tr n t:.~ 1 r 
O-Jr 2 - -J ■■ J 

300 

994.800 

994 .. 786 

493.. 791 

82. 068 

65.871 

/ c- '~4 r. ■ '**1 

oj n 

368 

1 230. 000 

1229. 963 

664. 203 

5 1 ! .. 34 1 

65. 935 

65 . 264 

420 

1478. 400 

1478.355 

827 . 058 

163. 404 

67 .477 

65. 292 

480 

1614.300 

1 6 1. 4 , 26 1 

1 046 - 753 

278. 142 

70. 3.22 

65, 3' ■/' 5 

540 

1637., 900 

! C>-~' / m . 2 . 

-i 3-rPi ■» r\ t 
I '7 i 1. O-.-' 

475 , 643 

75. 749 

0 u %-. h til V) 

600 

1654., 900 

1 654 .878 

13? 6. 532 

63 7.. 2j3£j 

85., 428 

66 „ 037 

660 

1663. 600 

1 663 .580 

1 ~:vP: . '-,94 

823 „ 638 

5,03. 338 

6 . 6 ., 806 

720 

1 647 . 800 

1 64“ 7 . 704 

13 43- 357 

948. 4 1 1 

1 34 „ 277 

66 . 466? 

780 

1647.400 

H / /I “T “T O /] 

i. Cj *+ / „ 00 ‘t 

6 361 E 6 3. 

3014.293 

176.813 

7 3. . 2 1 0 

848 

1652. 600 

.3 652.533 

1 378 , 793 

1 045 . 349 

O 0 0 O •“* 

75, 444 

988 

1 658 .400 

1 658 . 385 

3.39 5 . 9 46 

5. 05 7 . 3. 5 4 

y /j ”T “T 

\ , . A 

960 

1655.900 

1655, 837 

1 398 , 75 1 

5. 00,9 . 3 1 7 

297 ,. 873 

83 . 7, 60 

1 020 

1 477.600 

1477.606 

1 375 , 849 

1055.764 

323. 3. 15 

96 . 65 1 

1.080 

1293.000 

1 293 .012 

1. 299 . 525 

1036. 068 

339., 774 

3 05 . 760 

:l 1 40 

1 870.800 

1070. 821 

1 197*. 31 0 

995. 3 3 5 

348. 373 

1 1 5.339 

1 700 

9 I 5 . 4 88 

915. 41 8 

5 06? .. 564 

9?;5 .. 26 1 

347,975 

.! 25 . 032 

1 268 

809.700 

889 .. " 1 5 

•■'57. 899 

8.63 .. 837 

340 .. 490 

1 34 „ 47 ; 

1 328 

731. 988 

73 1 .. 9 1 2 

867 .. 8:54 

798 754 

328 , 799 

3 4 3 . 339 

1 3.80 

670 „ 708 

670. 71 0 

. 782.523 

72 3 . 88 8 

315. 724 

3 5 3. . 4 20 

1 44 0 

628 . 800 

620 888 

“ 1 . .. 288 

659 . 603 

383. 008 

3 53 „ 666 

1 580 

579 .. 788 

579. 787 

656 . 1 3 7 

604. 682 

29 3 . 44 7 . 

3 65. ■ 05 

1560 

544., 503 

544 ., 535 

c 84 „ 38 3 

[ j Cj » ShJL 5 ) 

28:1. ., 328 

1 70 .. 8 £5 

1620 

5 1 4 . 600 

5 1 4 684 

563 . 68 3 

515. 592 

272 . 66 / 

375,843 


Note: See figure 1 for locations of temperature simulations. 
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THEORETICAL MODELING 


The results obtained from both the heat chamber test and 
computer simulation show that the peak aluminum temperature of 
approximately 180 degrees F was well below the allowable maximum 
of 350. Since the aluminum temperature is the critical area of 
interest, it is desirable to know the optimum insulation thick- 
ness-peak aluminum temperature factor. To obtain this, ten 
additional computer simulations were conducted, each using a 
decreased total insulation thickness (rate of reduction was 
0.0095 feet per simulation). The results are presented in Figure 
8 as an insulation thickness versus peak aluminum temperature 
operating curve. 

This optimization run shows that the tile thickness could be 
reduced from 2.25 to 1.11 inches without exceeding the allowable 
350 degrees F structure temperature. Based on the materials used 
in the simulation, a 4 tile array which is 2.25 inches thick with 
a top surface area of one square foot weighs 1.1 pounds. The 
reduced insulation thickness tile weighs 0.65 pounds, resulting 
in a weight reduction of about 40%. 



PEAK ALUMINUM TEMPERATURE degree 



TOTAL INSULATION THICKNESS inches 


Figure 8. Theoretical Operating Curve 


CONCLUSIONS 


Within the framework of the governing assumptions, TEFS will 
provide a fast, low cost, accurate simulation of transient heat 
flow in a composite solid. As a aid in model design, TEFS 
enables a multitude of configurations to be analyzed in a short 
period of time, thus significantly reducing actual construction 
and testing costs. 
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